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ABSTRACT 

For half a century, evidence has been growing that the formation of stars follows a 
universal distribution of stellar masses. In fact, no stellar population has been found 
showing a systematic deviation from the canonical initial mass function (IMF) found 
for example for the stars in the solar neighbourhood. The only exception may be 
the young stellar discs in the Galactic Centre, which have been argued to exhibit a 
top-heavy IMF. 

Here we discuss the question whether the extreme circumstances in the centre of 
the Milky Way may be the reason for a significant variation of the IMF. By means of 
stellar evolution models using different codes, we show that the observed luminosity 
in the central parsec is too high to be explained by a long-standing top-heavy IMF as 
suggested by other authors, considering the limited amount of mass inferred from stel- 
lar kinematics in this region. In contrast, continuous star formation over the Galaxy's 
lifetime following a canonical IMF results in a mass-to-light ratio and a total mass of 
stellar black holes (SBHs) consistent with the observations. Furthermore, these SBHs 
migrate towards the centre due to dynamical friction, turning the cusp of visible stars 
into a core as observed in the Galactic Centre. For the first time here we explain the 
luminosity and dynamical mass of the central cluster and both the presence and extent 
of the observed core, since the number of SBHs expected from a canonical IMF is just 
enough to make up for the missing luminous mass. 

We conclude that observations of the Galactic Centre are well consistent with con- 
tinuous star formation following the canonical IMF and do not suggest a systematic 
variation as a result of the region's properties such as high density, metallicity, strong 
tidal field etc. If the young stellar discs prove to follow a top-heavy IMF, the circum- 
stances that led to their formation must be very rare, since these have not affected 
most of the central cluster. 

Key words: black hole physics - stars: formation - stars: luminosity function, mass 
function - Galaxy: centre. 



1 INTRODUCTION 



In ll955l . Salpeter found that the initial mass distribution of 
field stars in the range O.4<M*/M0<1O is a power-law with 
exponent 2.35. Since then, a large number of publications 
have investigated the initial mass function (IMF) of stars 
and made clear that star formation in general follow s the 
same empirical law, the canonical IMF |Kroupa 200l], and 
references therein). 

Due to its extreme conditions (mass density, veloc- 
ity dispersion, tidal forces), the Galactic Centre provides a 
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unique environment for testing the universality of the IMF. 
Star formation in the central region has thus been studied in 
detail, however no agreement has been reached on the nature 
of the IMF in either theory or observations: iManess et al.l 
(2007) find a best fit of observations in the central par- 
sec of our Galaxy with a model of constant star formation 
with a top-heavy IMF. On the other hand, iBuchholz et al.l 
(2009) show that the old stellar cluster in the Galactic Cen- 
tre very well resembles the bulge population. Observations 
of the young, massive Arches cluster in the central region 
of the Milky Way have long been interpre ted as a prime ex- 



ample for top-heav y star form ation ( e.g. Figer et al 
IStolte et all |2002|; iKim et al.l [200r3 : iKlessen et all 
However, lEspinoza et al.l (|2009T ) 



1999; 



20071 ). 



have shown that a canoni- 
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Table 1. Dependence of the composition of the Galactic Centre after 13 Gyr of star formation on IMF and star formation history. The 
columns give the following in order: Model name; underlying IMF (slope); range of stellar masses used; star formation history parameter; 
mass fraction in neutron stars; mass fraction in black holes; number of neutron stars and stellar black holes per 1.5 X 10 6 M (i.e. within 
1 pc from Sgr A*); number of bright stars (mag Ks < 17.5); mass-to-light ratio of unresolved population. While most stars survive the age 
of the Galaxy in the canonical model (a), an old cluster based on a top-heavy IMF is mass-dominated by stellar black holes. In contrast 
to the number of remnants, the mass-to-light ratio as well as the number of bright stars strongly depend on the star formation history. 



cal IM F cannot be excluded for this cluster. iPaumard et alj 
(2006) suggested a flat IMF for the young OB-stars observed 
in discs in the central parsec from the analysis of the K-band 
luminosity function. Based on m ore recent spectroscopic ob- 
servations, [Ba^c^^^^ J2009t]]find strong evidence for this 
to be true. iBonnell fc Rice! i|2008) found from SPH simula- 
tions that the IMF of stars forming in fragmenting accretion 
discs strongly depends on the parameters of the underlying 
gas infall scenario. Unfortunately, theoretical IMF predic- 
tions have failed in the past to co rrectly describe the obser- 
vations near the Galactic Centre (|Kroupall200"8ah . 

In this paper, we combine observational data with mod- 
els of stellar evolution and dynamics to constrain the stellar 
mass function and star formation history in the Galactic 
Centre. It is organised as follows: In Section [3] we analyse 
the properties of models of the Galactic Centre assuming 
different star formation histories, and compare them to the 
observations. Section [3] describes the mass profile of the 
central parsec and the effect of mass segregation. We dis- 
cuss the IMF of the young stellar discs around Sgr A* and 
appropriate formation scenarios in Section!?] and summarise 
in Section [5] 



2 STELLAR EVOLUTION MODELS OF THE 
GALACTIC CENTRE 

Observations of the central parsec of the Milky Way show 
that this region is dominated by a dense pop ulation of old 
stars with a total mass of ~ 1.5 x 10 6 M^ fGenzel et al.| 
20034 ISchodel et al. | 120071 , 120091 : iBuchholz et all l2009h. 
Within the uncertainties of a factor of two, Schodel et alj 



(2009) find that the extended mass inferred from kinemat- 
ics can be explained well by the visible stars. On the other 
hand, if star formation in the Galactic Centre occurs follow- 
ing a top-heavy IMF as suggested bv lManess et alj l|2007l ). 
one would expect a large number (and thus significant mass) 
of dark remnants. 

To test which mass functions are consistent with the 
observations, we u sed the stellar evolution package SSE 
|Hurlev et alJ koOO) to calculate population synthesis mod- 
els after 13 Gyr of star formation. Assuming (broken) power- 
law IMFs of the form £ (m) oc m~ a , where £(m)dm is the 
number of stars in the mass interval m to m + dm, we used 
the following models: 



(a) The canonical IMF according to iKroupal (|200ll ). 
£(m) oc m~ ai , with a = 0.3 (0.01 < m/M < 0.08), 
ai = 1.3 (0.08 < m/M < 0.5), and a 2 = 2.3 
(0.5 < m/M < 120). 

(b) A flat IMF with a = 1.35 (1 < m/M < 120) as sug- 
gested bv lPaumard et ail l|2006j ) for the young stellar discs 
in the Galactic Centre. 

(c) The same IMF, but exten ded to 0.01 < m/M < 120. 

(d) a — 0.85 as suggested bv lManess et alj ([2007) , again 
for 0.01 < m/M < 120. 

Starting with solar metallicity (Z = 0.02), we calculated 
evolutionary tracks for stars with masses 0.01 < m/M < 
120 in steps of 0.01 M . We averaged the results over time, 
weighted with different star formation histories of the form 
SFR(i) oc e - */* : We used constant (to = oo), exponentially 
declining (to = 3 Gyr, to = 1 Gyr, to = 300 Myr) and ex- 
ponentially increasing star formation rates (to = —3 Gyr). 
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Weighing the outcome with any of the IMFs above gives 
the total mass fraction of neutron stars (NSs) and stellar 
mass black holes (SBHs), as well as the total numbers of 
NSs and SBHs in the central parsec (assu ming an enclosed 
mass of 1.5 x 10 6 Mm : ISchodel et alj|2009h . To estimate the 
total number of bright stars (mag Ks < 17.5) and the K-band 
mass-to-light ratio M / Lk b of the unresolved stars and stel- 
lar remnants, we generated a sequence o f Padova isochrones 
ijMarigo et al.ll2008l : iBertelli et al.lll994l ). assuming an aver- 
age extinction of 3. 3 mag for stars in the Galactic Centre 
l|Schodel et al.l [2009). We used the magnitu des calculated 
for the 2MASS Ks filter ijCohen et al.ll2003h . whose trans- 
mission curve closely resembles that of the Ks-band filter of 
the NAOS-CONICA instrument at the ESO VLT, which was 
used for the observations of the Galactic Centre discussed 
here. We further assumed that the mass of a stellar remnant 
depends on the initial mass as 

f 0.109 mjnit + 0.394 M , 0.8 < m init /M < 8, 
m Tem = I 1.35 M Q 8 < m in i t /M Q < 25, 

(o.lm init , 25<m init /M . 



dKalirai et all 20081: Tho 


rsett & Chakrabartvl 1999 


iBaumeardt & Mieskd 120081; 


Dabringhauscn et al. 2009) 



The uncertainties of these assumptions are not easily 
quantified, since the underlying theory and observations 
are not robust, but using the masses, radii, and effective 
temperatures from our SSE models and assuming a black- 
body spectrum leads to comparable results. We are thus 
confident that our results are correct within a factor much 
less than two. 

Table [T] lists the respective values for all four models. 
While most stars survive the age of the Galaxy as main 
sequence stars in the canonical model (a), an old clus- 
ter based on a t op-heavy IMF is ma ss-dominated by stel- 
lar black holes. ISchodel et ail {2009) find that the unre- 
solved stellar population makes up > 98% of the mass in 
the central parsec, and find its mass-to-light ratio to be 
M/L Ks = lAtg* M /L 0iKs - Figure [T] plots the ratio of to- 
tal mass to diffuse light in our models as a function of the 
star formation history parameter to. We find that the ob- 
servations are consistent with the canonical IMF (a), with 
a tendency towards constant or increasing star formation. 
The a = 1.35 models (b, c) require increasing star forma- 
tion, and an IMF as flat as a = 0.85 (d) is not consistent 
with the observed old population at all. 

We can also compare the results of our analysis as given 
in Table[T]to stellar number counts in the central parsec. As- 
suming a canonical IMF, we expect ~ 2.5 x 10 4 SBHs and 
NSs for every 1.5 x 10 6 M in stars and stellar remnants, the 
latter being the estimated enclosed mass within 1 pc from 
Sgr A*. Due to dynamical friction, SBHs may migrate to 
the central parsec from as far as 6 pc from the centre within 
10 Gyr (see Section [3] for details), thus increasing the num- 
ber of SBHs in the central parsec by up to one order of 
magnitude. On the other hand, SBHs may spiral into the 
SMBH as discussed above . In total, t h e exp ected number 
of SBHs agrees best with iMuno et alj (|2005l ) suggesting a 
number of ~ 10 4 SBHs and NSs in this region from X-ray 
observations, while models based on a top-heavy IMF sug- 
gest numbers well above 10 5 . 

A more verifiable quantity is the number of bright stars: 
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Figure 1. Ratio of total mass and diffuse light in our models as a 
function of the star formation history para meter tp. The s haded 
area marks the one-sigma range derived bv lSchodel et all (120091 ) 
from observations. While observations are best explained with an 
increasing star formation rate following the canonical (a) or a 
moderately top-heavy (b, c; a = 1.35) IMF, observations cannot 
be reproduced assuming a flatter IMF (d; a = 0.85) 



ISchodel et ail (|2009h find ~ 6000 stars with Ks-band mag- 
nitude <17.5 within a projected distance from Sgr A* of 
1 pc. Assuming a spherical distribu tion with a density pro- 
file p X r" 175 (|Schodel et al.ll2007D . ~ 3000 of these stars 
are in the innermost parsec. This value is consistent with our 
model (a) of a canonical IMF, assuming a constant or declin- 
ing star formation rate, or with models (b) and (c) assuming 
constant or slightly decreasing star formation (to > 3 Gyr). 
An even flatter IMF (d) requires a star formation rate in- 
creasing with time to explain the observed number of lumi- 
nou s stars. 

iBuchholz etafl (2009) find that the K-band luminos- 
ity function (LF) of late-type stars in the central parsec 
closely resembles that of the bulge population. Here we use 
the luminosities calculated from our stellar evolution models 
to compare the LFs of different IMFs. Figure [2] shows the 
K-band LFs of our models as a function of star formation 
history. It is seen that the shape of the LF does not depend 
significantly on the IMF. In particular, all curves can be ap- 
proximated by a power-law with a slope j3 ~ 0.3, and all 
exhibit the horizontal branch / red clump peak. Th is peak 
is offset to the observations of lBuchholz et all (|2009h by one 
magnitude, which may be due to our models of stellar evolu- 
tion and the assumed constant extinction. We thus presume 
an uncertainty of a factor of two in the luminosity estimates 
discussed above, which does not affect our qualitative results 
and reasoning. 

Altogether, observations of the old stars, diffuse light, 
and stellar dynamics in the Galactic Centre are best ex- 
plained with a canonical IMF at constant or decreasing star 
formation, but may also be explained with a somewhat flat- 
ter IMF and an almost constant star formation rate. An 
IMF slope flatter than a = 1 can be safely ruled out. 

We have to stress that here (and in the next section) 
we only calculate self-contained models, i.e. we assume that 
the stars now present in the Galactic Centre also formed 
in this region. Stars may also be brought to the vicinity of 
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Figure 2. Lum inosity function of brig ht stars in the Galactic Centre, assuming a K-band extinction of 3.3 mag. The panels compare the 
LF observed by Buchholz et al. (2009|), represented by their fit as a straight line (and corrected for projection effects), to the results of 
our models of different IMFs following exponentially increasing (upper left), constant (upper right), and exponentially decreasing star 
formation (lower panels). The slope derived from the observations fits all our models quite well: It does not depend strongly on star 
formation history or IMF. However, as is seen also from Table [JJ the normalisation does depend on the star formation rate, especially so 
for the top-heavy models (d, b). 



the SMBH by capture of individual stars or clusters, and 
some unknown fraction of the stars observed in the centre 
of the Milky Way may have formed in a much different en- 
vironment. To this extent, we do not exactly discuss star 
formation in the Galactic Centre, but the formation of stars 
now observed in the central region, thus reaching an insight 
into the composition of stellar systems surrounding SMBHs. 



3 MASS PROFILE IN THE GALACTIC 
CENTRE 

Theoretical arguments and TV-body simulations show 
that a stellar system around a SMBH evolves into 
a cusp with a 7 = 1.75 power-law density dist ri- 
bution (iBahcall fc Wolj 1 19761; IBaumgardt et ail l2004al lbt 
IPreto et alj 120041 ; lAmaro-Seoane et all |2004 ). Indeed, ob- 
servations show that the central pars ec of the Milky Way 
exhibits a corresponding profile (e.g. iGenzel et al. 2003b; 
ISchodel et alj [20071 . l2009h . However, observations also re- 
veal a deficit of stars within a few 0.1 pc fr o m the SMBH 
llGenzel et al.lll996l . [2003al ; 
120071. 



2009): For example, 



Figer et al J 120031; ISchodel et all 



Schodel et al.1 (|2007f ) estimated 



the density profile of the stellar cusp in the Galactic Centre 
from the observed luminosity profile and kinematics as 



p(r) = (2.8 ± 1.3) x 10 6 M pc" 



0.22 pc 



where r is the distance from the SMBH, 



7 



(2) 



1.2 inside 



0.22 pc and 7 = 1.75 outside 0.22 pc. iBuchholz et al.1 (|2009T l 
find that the profile may be flatter or even slightly inverted 
in the central 0.2 pc, turning the cusp into a core. 



The observed core in the old population around 
Sgr A* may be explaine d by stellar collisions destroying 



the envelopes of giant s ([Ge nzcl ct al. 1996; Davic s et al 
1 19981 ; lAlexanderl 1 19991 ; iBailev fc Daviesl Il999l ; iDale et al 
120091 ) . or by an inspirall i ng IMBH depleting the central 
region ([Ba umgard t et ail 120061 ; lLockmann fc Baumgardt] 
l2008l ; lMerritt fc Szellll2006l ). 

Here we suggest that the central 0.2 pc around Sgr A* 
are not mass-depleted, but dominated by a cusp of SBHs 
having displaced the less massive visible stars: Due to dy- 
namical friction, the most massive stellar remnants sink into 
the depth of the SMBH's potential well. Thus, the SMBH 
is expected to be surrounded by a large number of black 
holes in its imme di ate vicinity (IBaumgardt et al.l [2004b; 
iFreitag et al] 120061 ). iFreitag et al.l (|2006f t have shown that 
the SBHs in the Galactic Centre establish a cusp profile 
p(r) oc r -1,75 , while the density profile of the less massive 
stars becomes consistent with the observed luminosity pro- 
file. They do not find a clear density cutoff around 0.2 pc as 
in the observations, since they had to use an unrealistically 
large number of SBHs due to computational limitations of 
the methods used. Since it is as yet not possible to perform 
a full calculation including a realistic number of stars and 
SBHs, we here follow a theoretical approach. 

For simplicity, we assume that the Galactic Centre 
formed evolved stars and SBHs 10 Gyr ago without any fur- 
ther star formation. The frictional drag on an inspiralling 
SBH can be estimated as 



dv 
dt" 



4tt In AG 2 p*(r,t)M B 



n V 

erf(X) - =^e 
V 71 " 



(3) 



l|Binnev fc Tremaind[l987l . Eq. 7-18), where v is the veloc- 
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ity of the SBH, G is the gravitational constant, AIbh is the 
SBH mass, p*(r, t) is the density of the stellar background, 
In A is the Coulomb logarithm and X = v/(y2a) is the ra- 
tio between the SBH velocity and the (ID) stellar velocity 
dispersion a. For a stellar density profile p*(r) oc r ' with 
1-2 < 7 < 1.75, one gets X ~ 1.2. Assuming that the over- 
all density profile does not change with time and can be 
described by p(r) = p*(r,t) + pBTl(r,t) = por 1 yields 



s (r) = <G 



MSMBH + 7! r 

3-7 



3-7 



(4) 



as the circular velocity at distance r of the SMBH of mass 
Msmbh. Inserting into 



dv 
dt 



Ft 
Mbh 



(•») 



where F and L are the frictional force and orbital angular 
momentum, respectively, leads to the SBH inspiral speed 



-87rr 5 / 2 lnAv/Gp*(r,t)M B H [erf(X) - f§e~ x2 ] 
( Msmbh +4np ^r 3 -^ ^Msmbh + f^r 3 ""? 



(6) 



assuming a circular orbit of the SBH. Note that p*(r, t) de- 
creases with time due to the inspiralling SBHs increasing 
PbhM). 

To integrate this equation numerically, we assume that 
the cluster was initially not mass-segregated, pbh(j*, 0) = 
0.04p(r), where Mmi/Mtot ~ 0.04 was taken from 

Table [TJ for a canonical I MF, and p(r) = 2.8 x 

lO 6 M pc- 3 (r/0.22pc)" 7 is a iBahcall fc Wolj l|l97rj | pro- 
file with 7 = 1.75. We further assume a SMBH mass of 
Msmbh = 4 x 1O 6 M , and SBHs of mass Mbh = 10 Mq. 
Figure [3] shows the enclosed mass in SBHs as a function of 
central distance, Mg B H(< r ), for different times t in steps of 
1 Gyr, starting with Mg B H(< r ) oc r 3-7 at t = 0. It is seen 
that the cusp is saturated by SBHs in the innermost 0.5 pc 
within 10 Gyr. Clearly, this is not the final answer: Due to 
random deflections of stars, we can expect a number of stars 
in the innermost region. Furthermore, contrary to our as- 
sumptions above, both stars and SBHs move on eccentric 
orbits, preve nting a strict s egreg ation. 

In fact, Frei tag et all (120061) find that the inspiralling 



SBHs build a lBahcall fc Wolj (|1979 ) profile, while the den- 



sity profile of the stars is flatter, compatible with the slope 
observed in the central 0.2 pc. In particular, they find that 
the central part will be mass-dominated by the SBHs, sug- 
gesting that the estimates above are at least qualitatively 
correct. 

Altogether, the above discussion suggests that 

• the innermost region is mass-dominated by SBHs within 
a ~ 0.5 pc, the semi-major axes a being distributed as 
a 2 " 1 - 75 , 

• in this region, the stars are distributed roughly as 
a 2-1 ' 2 , and 

• outside a ~ 0.5 pc, the cluster is strongly dominated by 
visible stars, p* > 0.94p assuming a canonical IMF. 

To compare these results to the density profile derived from 
observations, we need to convert the distributions of semi- 
major axes into a distance distribution. For this, we assume a 
thermal distribution (i.e. a uniform distribution in e 2 , where 
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Figure 3. Mass segregation in the Galactic Centre. The solid 
lines show the enclosed mass in SBHs as a function of central 
distance in steps of 1 Gyr from the black bottom curve (f = 0) 
to the bright top curve (t = 13 Gyr). The dashed line depicts the 
total enclosed mass in SBHs and stars, which is assumed to be 
constant. It can be seen that the central 0.5 pc are saturated with 
SBHs over the Galaxy's lifetime. 



the eccentricity e is distributed as / e (e) = 2e). A cusp with 
3D "density" of semi-major axes p a {a) has a ID distribution 
4>a(a) = A-Ka 2 p a {a). The central distance of a particle on a 
Kepler orbit with semi-major axis a and eccentricity e passes 
any value r (with 1 — e<r/a<l + e) twice during a full 
orbital period T, hence its distribution function is 



fr{r) 



{a-rf 



(7) 



A particle on an orbit of eccentricity e and semi-major axis 
a assumes central distances r between pericentre and apoc- 
entre distance, a(l — e) < r < a(l + e). Depending on the 
orbital eccentricity, a particle at distance r can thus have a 
semi-major axis between r/2 and infinity. For a given semi- 
major axis a, the eccentricity has to be larger than |1 — r/a\ 
to be consistent with the central distance r. Altogether, the 
distribution function of central distances is 



<t>r{r) 



f e (e)cf)a(a)f r (r) de da 



a — r/2 J e— 11 — r/a 
1 



8a r ep a (a) 



; de da, 



l a-r/U e=\l-r/a\ ^a 2 e 2 — (a — r) 2 
and with (f> r (r) — 4ivr 2 p(r), the 3D density profile is 

2a e p a (a) 



p(r) = 



o=r/2 Je=\l-r/a\ 7rr^/a 2 e 2 — (a — r) 2 



; de da. 



(8) 
(9) 

(10) 



Figure U shows a semi-major axis distribution p*, a (a) = 
2.8 x 10 6 M Q pc~ 3 (0.5/0.22)- 7 x 0.96 of stars in the Galac- 
tic Centre, as well as the resulting 3D density profile p*(r). 
It can be seen that the density profile derived in our the- 
ory resembles the observations very well. The exact value 
of the break radius cannot be determined within the uncer- 
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x I pc 



Figure 4. Density and semi-major axis distributions of stars in 
the Galactic Centre. The dashed line shows a semi-major axis 
distribution p*, a (a) with a break radius of 0.5pc, as predicted by 
theory for a canonical IMF. Assuming a thermal eccentricity dis- 
tribution, the dash-dotted line shows the corresponding density 
profile p *(r), which resemble s the profile derived from observa- 
tions by Scho del et al ] l|2007l l very well. Both distributions are 
consistent with the observed profile within the lc uncertainty 
shown as a shaded area. In contrast, the expected density pro- 
file for a top-heavy IMF (model (c) in Table [TJ dotted line) is 
not consistent with the observations, producing a break radius at 
2.5 pc already after 8 Gyr of mass segregation. 



tai nties of the density p rofile derived from the observations 



tamties ot tne den sity p : 
bv lSchodel et~ai1 l|20071 ). 



On the other hand, models with a black hole mass ratio 
MBH/Aftot ^ 0.6, as would result from the evolution of a 
population of stars following a top-heavy IMF (cf. Table [TJ, 
lead to a break radius of 2.5 pc and higher already after 
8 Gyr of mass segregation. Thus, the observed break radius 
is further evidence for a canonical IMF and a moderate black 
hole fraction in the Galactic Centre. 



4 IMF OF THE YOUNG STELLAR DISCS 

Observations of the Galactic Centre revealed one or 
two discs of ~ 6 Myr old stars orbiting the central 
super- massive black hole (SMBH) at a distance of ~ 
0.1 pc llLevin fc Beloborodovl 120031 iPauniard et al.l l200fj ; 
iLu et al.l l200d. 120091; iBartko et all l2009al ; see discussion 
in lLockmann fc Baumgardt! |2009| ). Two main scenarios 
have been proposed for the formation of such massive 
(Mdisc ~ 1O 4 M ) discs close to Sgr A*: iGerhardl (|200ll ) 
suggested that a disc of stars may have formed by tidal 
disruption of an infalling cluster of young stars, which 
would require a large mass or a central intermediate- 
mass black hole to survive the strong tidal forces from 
the SMBH (iMcMillan fc Portegies ZwartJ I2003T ) . However, 
iLevin fc Beloborodovl (|2003l ) showed that this cannot ex- 
plain the small distances of the stars from Sgr A*, 
and propose in-situ formation by fragmentation of a 
mas sive accretion disc as an alternative sc e nario (see 
also INavakshin fc Cuadral 120051 ). iMapelli et all l|2008h and 
iBonnell fc Rice! ( 20081 ) have shown the infall of a giant 
molecular cloud (GMC) towards the Galactic Centre to 



be effective in creating a disc of stars with distances from 
Sgr A* consistent w i th th e observed young massive stars. 

iPaumard et al.l (|2006l ) derive a flat IMF for the observed 
disc stars from their K-band luminosity function; however, 
this only refers to a mass interval of 20 < M+/Mq < 30. 
From th e amount of mutual wa rping of the observed stel- 
lar discs, INavakshin et al] (|2006l ) derived upper mass limits 
to the discs' masses consider i ng th eir apparent flatness as 
measured by IPaumard et al. (120061) after a f ew Myr of in- 
teraction. However, in lLockmann et al] |2009l) we show that 
a canonical IMF cann ot be excluded from disc dynamics: 
INavakshin et al] £2006) used stellar discs only out to 0.2 pc 
fro m the SMBH, while mo re than one third of the stars listed 
by IPaumard et al. I (|2006l ) are further away from Sgr A*. In 
addition, they did not distinguish between stars in the outer 
parts of the discs which almost retain their orbital planes, 
and stars close to the centre whic h strongly precess. 

INavakshin fc Sunvaevl l|2005h argue that the X-ray lu- 
minosity of the Sgr A* field is too low to account for the 
number of young < 3 Mq stars expected from a canonical 
IMF, considering the large number of O-stars observed in 
the discs, which may be explained by a higher low-mass 
cutoff near 1 Mq . On the other hand, the existence of the 
S-stars within 0.01 pc from Sgr A* is in favour of a canonical 
IMF (and thus a large number of B-type stars) for the stel- 
lar dlsc^ J Jftiisy_wreJomied from these discs as suggested 
by lLockmann et al] (|2008h . For a more de tailed discussion 
of the IMF of the young stellar discs, see lLockmann et al] 
(|2009l ). 

Assuming an infalling cluster as the origin for the stellar 
discs, one may expect a top-heavy mass function if the clus- 
ter was mass segregated and then tidally stripped. On the 
other hand, a top-heavy initial mass function of stars formed 
in a fragmenting disc would require an unusual mode of star 
formation. As we have discussed in the previous sections, no 
convincing reason or evidence for a flat IMF in the Galactic 
Centre in general has been found. 

Using SPH si mulations of star f orma tion in fragmenting 
gas accretion discs, IBonnell fc Ricel (|2008T l find that the IMF 
of disc stars can be bimodal (and thus top-heavy) if the 
infalling gas cloud is massive enough (>10 5 M Q ) and the 
impact parameter of the SMBH encounter is as small as 
~ 0.1 pc. This way, an extreme configuration of cloud masses 
and distances may lead to a significant variation of the IMF. 

Only recently, a systematic search of OB stars in the 
central parsec revealed a significant deficit of B-type stars 
in the regime of the young discs, suggesting a strongly top- 
heavy IMF for these discs (|Bartko et al]|2009bl ). Until then, 
there was no strong evidence for a top-heavy IMF in the 
young discs. However, it cannot be assumed that the ma- 
jority of stars in the Galactic Centre were formed the same 
way as the young disc stars, whose existence may be an in- 
dication of recently enhanced star formation processes: If 
the Galactic Bar is young, as an incr eased fraction of b arred 
galaxies for lower redshifts suggests (|Sheth et alj|2008l ). bar- 
induced gas inflow may explain such an enhancement by an 
increas ing supply of high-mass GM Cs towards the central 
region (|Sellwood fc W ilkinson 1993! and references therein). 
Hence, despite the significant lack of B-stars in the range 
0.03 — 0.5 pc from the SMBH, suggesting that the young 
discs indeed formed following a flat IMF, this does not im- 
ply that star formation in the Galactic Centre is or was in 
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general top-heavy. Instead, the mass function of disc stars 
may reveal deta i ls of t he formation scenario, as the results of 
iBonnell fc Rice! HH) suggest. The majority of old stars in 
the region may thus have formed under different conditions 
(e.g. mass and impact parameter of the infalling clouds) in 
the scenario of a fragmenting disc, or formed further away 
and then migrated to the centre (as in the infalling cluster 
scenario, see above), or formed by any other process follow- 
ing the canonical IMF. 



5 DISCUSSION 

Various attempts have been made to study star formation in 
the Galactic Centre, but so far neither theory or simulations 
nor observations led to an agreement on the (initial) distri- 
bution of stellar masses. Here we have shown that theory 
and observations are consistent with star for mation gener- 
ally following a canonical IMF (|Kroupall200ll ') in the Galac- 
tic Centre, just as anywhere else in the Universe. Our main 
results can be summarised as follows: 

1. The mass-to- light ratio of the central parsec of the 
Milky Way is consistent with a constant or exponentially 
decreasing star formation rate following a canonical IMF. 
Models of constant star formation following an IMF with 
q = 1.35 are consistent with the observed luminosities but 
create ~ 10 5 SBHs in the central parsec, ten times more 
than expected by other authors. Mass functions flatter than 
a ~ 1 can be safely ruled out, since they cannot explain the 
observed number of bright stars and the diffuse light. 

2. The core observed in the luminosity distribution with 
a radius of rb rC ak ~ 0.2 pc does not imply a core in the 
mass profile, bu t can be well explain ed by mass segregation 
as suggested bv lFreitag et al.l (|2006l ). where dark remnants 
mass-dominate this region. Again, the observations are best 
explained by a canonical IMF, and are not compatible with 
star formation following a top-heavy IMF with a < 1.35, 
as this would create a core radius one order of magnitude 
larger. 

3. Recent observations revealing a deficit of B-type stars 
in the young stellar discs suggest a top-heavy IMF for this 
population, which may be explained by tidal stripping of an 
infalling mass-segregated cluster, or unusual modes of star 
formation in a fragmenting accretion disc. However, these 
results do not allow conclusions on star formation in the 
Galactic Centre in general, for which we have no reason to 
assume it to be non-canonical. 

While other au thors generally predicted a flat I MF (e.g. 
Klessen et al 1I2007TI or a higher lo w- mass cut-of f (e.g. lMorrisl 
19931 ; iLevin fc Beloborod^vll2003l ; lLa7sonll2006l ) from state- 
of-the-art theoretical star formation models of the Galactic 
Centre, we find that observations suggest star formation fol- 
lows a canonical IMF even under the extreme circumstances 
present in the central cluster. This universality of the IMF 
poses a major challenge to our understandin g of star forma- 
tion processes (see also lKroupall200ll . 12008b) . 

Our results rely on the assumption that the stars ob- 
served within 1 pc from the Galactic Centre also formed 
there, suggesting star formation with a canonical IMF even 
under such exotic conditions. It is possible that some of 
the stars were brought in by massive star clusters which 



spiralled towards the SMBH through dynam ical friction 
jPortegies Zwart et all 120061 ; iFuiii et aljfeoOSl) . However, 
this scenario is unlikely because a star cluster is stripped 
on its way towards the centre and loses mostly low-mass 
stars, since it would be in a mass segregated state soon after 
its formation. Therefore, the most likely scenario is a cen- 
tral cluster that formed over a Hubble time with a canonical 
IMF, where the very young stellar population of the stellar 
discs o bserved to have a very top-heavy IMF (|Bartko et alj 
2009b) constitutes a rare star formation event not typical 
for the bulk stellar population in the central cluster. 
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